Correlation of transverse momentum and multiplicity 
in a superposition model of nucleus— nucleus collisions 
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In p-p collisions the average transverse momentum is known to be correlated with the multiplicity 
of produced particles. The correlation is shown to survive in a superposition model of nucleus- 
nucleus collisions. When properly parameterized, the correlation strength appears to be independent 
of the collision centrality - it is the same in p-p and central A-A collisions. However, the correlation 
is strongly suppressed by the centrality fluctuations. 
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(NJ , I. INTRODUCTION 

The transverse momentum of particles produced in p-p collisions is correlated with the particle's multiplicity. At 
' the SPS energy the correlation, which is directly observed experimentally 00, is most probably of simple kinematical 
origin - when the multiplicity N of produced particles at fixed collision energy grows, there is less and less energy 
' to be distributed among transverse degrees of freedom of produced particles. Consequently, the average transverse 
momentum {pt)(n) decreases when N grows. The dependence of {pt)(n) on N is approximately linear QLQ- 

The correlation px vs. N seems to dominate the transverse momentum fluctuations which have been experimentally 
studied [3] in nucleus-nucleus (A-A) collisions at SPS energy. The fluctuations have been quantified by the measure 
$(pt) 0] and it has been shown 0] that the non-monotonic behavior of $>(pr) as a function of collision centrality can 
be explained by the observed non-monotonic multiplicity fluctuations 0,0 combined with the correlation pt vs. N. It 
contradicts an expectation that the correlation, which is well seen in p-p interactions, is washed out in A-A collisions. 
To clarify the situation the NA49 Collaboration has undertaken an effort to directly measure the correlation pr vs. N 
\Q , in A-A collisions at precisely determined number of projectile participants. 

An analysis of the experimental data, which is underway, is performed for p-p and Pb-Pb collisions at 158 AGeV in 
the forward rapidity window [4.0,5.5] (the center-of-mass rapidity is 2.92). The collision centrality is determined on 
■ the event-by-event basis, using the zero-degree calorimeter which measures the energy deposited by remnants of the 
1— h ' projectile nucleus. The energy is recalculated into the number of spectator nucleons coming from the projectile. The 
experimental method was worked out in the course of multiplicity fluctuation measurements performed at precisely 
fixed collision centralities 0, . 

The motivation of our analysis is twofold. The first goal, which is oriented towards experimentalist's needs, is 
to check how a finite calorimeter resolution, which hinders the centrality determination, influences the observed 
correlation. The second goal is to formulate predictions within simple models of nucleus nucleus collisions as a 
reference for the experimental data. For these purposes we discuss the correlation px vs. N in A-A collisions, assuming 
that the particles produced in such a collision originate from several identical sources which are independent from 
each other. A particular realization is the model where nucleus-nucleus collision is a superposition of independent 
nucleon-nucleon interactions. Then, a single nucleon-nucleon interaction is a single source of particles. More realistic 
is the Wounded Nucleon Model Q where the wounded nucleons - those which interact at least once - contribute 
independently from each other to the final multiplicity. Then, there are two wounded nucleons in a nucleon-nucleon 
interaction, and thus, there are two particle's sources. Since the particles in the forward domain are analyzed by the 
NA49 Collaboration, only the wounded nucleons form a projectile are expected to matter, and then their number 
equals the number of particle's sources. 

Our results obtained within the approach of independent sources can be interpreted, referring to one of three 
categories of models of nucleus- nucleus collisions introduced in There are transparency models when the projectile 
(target) participants, which are excited in the course of interaction, mostly contribute to the particle production in 
the forward (backward) hemisphere. One deals with (rather unrealistic) reflection models when the projectile (target) 
participants mostly contribute to the backward (forward) hemisphere. Finally, there are mixing models when the 
projectile and target participants equally contribute to the backward and forward hemispheres. The latter models 
were shown to properly describe (large) multiplicity fluctuations of particles produced in the forward hemisphere 
at fixed number of projectile participants 00, as the target participants of fluctuating number contribute to the 
forward hemisphere. The transparency models highly underestimate the multiplicity fluctuations while the reflection 
models strongly overestimate them. 
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Because of various possible interpretations of our results, we use in Sees. IIIIIVI the term source not referring to 
a specific model. We first discuss the correlation p? vs. N for particles coming from a single source (Sec. [nj, and 
then, we consider particles which originate from several identical and independent sources fSec. IIII(I . In Sec. IIVI a 
situation with varying number of particle sources is discussed. Our results are interpreted in two models in Sec. Ivl 
where specific predictions are presented. To simplify the notation, the transverse momentum is denoted from now on 
as p not as pr or p±. The word 'transverse' is also mostly skipped in the text. 



II. SINGLE SOURCE 



Pn(P1iP2i ■■ -iPn) is the probability density that N particles with momenta pi, p%, . . .p^ are produced by a single 
source. It is normalized as 

^2 / d Pldp 2 ■ ■ ■ d PNPN(Pl,P2, ■ ■ ■ ,Pn) = 1 • (!) 
The multiplicity distribution Vn is given by 

Vn =' J dp 1 dp 2 ■ ■ ■ dp N P N (pi,p 2 , ■ ■ ■ ,Pn) ■ 

The momentum distribution of a single particle in the events of the multiplicity N is defined as 
def 1 1 



P ^ {p) p n n 



J dp 2 ■ . .dp N P {N) {p,p 2 , ■ ■ ■ ,Pn) + J dpidp 3 . ..dp N P( N ){pi,p,p 3 ■ ■ ■ ,Pn) (2) 
+ h J dpidp 2 ■ ■ .dpN-iP{N)(piiP2 ■ ■ ■ ,PN-1,P) , 



and it is also normalized 

dpP(N) ip) = 1 • 

The single particle average momentum in the events with N particles, which is denoted as (p)(n), is correlated 
with TV, and this correlation has to be parameterized for our further considerations. According to the data 0,0 > the 
dependence of (p)n\n on N is approximately linear, and thus we write 

(P>W = « + 6 ( 1 -^) ■ ( 3 ) 
We note that the parameter a only approximately equals the inclusive average momentum (p) which is 

<P> = (j^E^OV)' ( 4 ) 



N 



with 



The parameterization @ gives 



N 



(p) = a + 6- 



(N) 2 ' 

which for the case of poissonian multiplicity distribution is 

(p)=fl -(A7) ' 

The multiplicity distribution of charged particles produced in high-energy proton-proton collisions is known not 
to be poissonian. First of all, the total number of charged particles is always even due to the charge conservation. 
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The multiplicity distribution of positive or negative particles is not poissonian as well - the variance does not grow 
linearly with (N) but it follows the so-called Wroblewski formula *- e - the variance is a quadratic function of (N). 
However, for the collision energies of the SPS domain, which are of interest here, the Poisson distribution provides a 
reasonable approximation. The multiplicity distribution of not only positive or negative but of all charge particles is 
poissonian with high accuracy when the particles are registered in a limited acceptance. Therefore, the multiplicity 
distribution, which is further used in our calculation, is assumed to be of the Poisson form. 



III. SEVERAL SOURCES 



The momentum distribution of N particles coming from k identical and independent sources equals 

JVi N 2 N k 

x Pn 1 (pi,P2, ■ ■ • ,PnJ Pn 2 {pn 1 +i,Pn 1 +2, ■ ■ ■ ,Pn 1 +n 2 ) 

X • ■ • X P Nk (p Nl + ,,, + Nk _ 1 + i 1 p Nl + ,,, +Nk l+2 , ■ ■ ■ ,PNi+...+N h ) ■ 

One checks that the distribution (JjjJ is normalized in agreement with Eq. Q. One defines the single particle distri- 
bution analogous to and computes the average single particle momentum at fixed multiplicity as 

<P> W = W7f ££••■£ 5 n 1+N2+ - +N '"Pn 1 V N2 ■■■V Nk (Ni(p) w + N 2 (p) {N2} + ■■■ + N k (p) {Nk) ) . (6) 

N N± N 2 N k 

Using the parameterization J3J|, one finds 

<p>w - ^r^£E-£^ 1+%+ - + ^^-^ ( 7 ) 

N Ni N 2 N k 



/ / N 2 N 2 N 2 \ 

x ( W+ * + ... + * 4 )(. + t )- l (^ + ^ + ... + ^) 



V/ I N Nl N2 Nk 



k) > 



where it has been observed that 

Ni + N 2 + ■ ■ ■ + N k = N , 
(Ni) + (N 2 ) + ■ ■ ■ + (N k ) = (N) k , 

n = ££••• £ 4 M+ - +Nk ^VN 2 ■ ■ ■ V Nk . 

Ni N 2 N k 

The average with the single-source distribution is denoted, as previously, as (• • •) while the average with the k— source 
distribution is (■ ■ -} k . 

Further calculation is performed with the poissonian multiplicity distribution which is 

•-«•>■ 

The convolution of k identical Poisson distributions of the average (N) is known to be equal to the Poisson distribution 
of the average (N) k =k{N). Thus, 

££•■■£ S^ + - +N ^V Nl P N2 ...V Nk = e-<">» = V% . (8) 

Ni N 2 N k 

Using the same technique, which provides Eq. (JSJ), one proves that 

££•••£ W +N2+ - +Nk r Nl v N2 ...v Nk Nf = (f + ■ 0) 

Ni N 2 N k 
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Substituting the formula into Eq. Q, one finds the result 

W M =fl+i ( 1 "T)' (10) 

which gives Eq. (J3J) for fc = 1. If (TV) > 1 or equivalently (N)k 3> fc, the explicit dependence of (p)?jv) on & i s 
effectively very weak. For the Poisson distribution with (N)k ^> fc > 1, the multiplicities Af of sizable probabilities 
do not much differ than (N)k. Consequently, N + k — 1 can be approximated by N in Eq. and we get the result 
analogous to the k = 1 case (01 that is 

(^)-^(l-^). (ID 

Thus, the average momentum at fixed multiplicity is approximately independent of the number of sources as it depends 
only on the ratio N/(N)k- The correlation strength parameter b is entirely independent of the source number. 



IV. VARYING NUMBER OF SOURCES 



In the previous section the number of particle's sources has been treated as a fixed number. In a real experimental 
situation, however, the number of sources cannot be fully controlled in nucleus-nucleus collisions. The centrality 
selection is never perfect but even an exact determination of the participant number does not guarantee that the 
number of particle sources is always the same for a given number of participating nucleons. Thus, let us consider 
what happens when the number fc varies event by event while the multiplicity of produced particles is kept constant. 

When the number of particle sources fluctuates, the momentum distribution of N particles coming from fc sources, 
which is given by Eq. (JSJ, has to be averaged with the distribution of number of sources Pfc. The analogue of Eq. © 
reads 

«p>fo> = 7ir^E p fe EE---E< 1+Ar2+ - +Arfc ^^---^ d 2 ) 

^ N > k Ni N 2 N k 

x (Ni{p)( Nl ) + N 2 (p)(n 2 ) + ■■■ + N k (p)(N h )) , 
where (7^) is the multiplicity distribution averaged over number of sources which equals 

(K) = T,^T,T,---i: C +N2+ - +Nk -P Nl r N2 ...v Nk . 

k N\ N 2 N k 

Substituting the parameterization © into Eq. I)12[l . one finds 

«p>fo> = 7iy^E^EE---E^ 1+JV2+ - +JVfc ^^---^ 

\ n/ k Ni N 2 N k 

x {^^^■■■^ k){ a + b)-b{^ )+ ^ )+ ... + ^)) 

= - + * - 7^4^ E p**EE • • • E 4 1+N > + - +N ^ Nl v N2 ...v Nk ni . 

When the single-source multiplicity distribution Vn is poissonian, one uses the formula © and Eq. (|13|l gives 



(13) 



«*>(»)> = a + b ( l - T^a - 7^ E p 4 K) • (W) 



Analytic calculations can be easily performed when the number of sources is treated not as an integer but as a 
continuous (positive) number distributed according to the Gamma distribution 

P(*) = Y^/e^ ■ W 
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where T(x) is the gamma function and the parameters a, j3 obey a > and j3 > — 1. One easily finds the average 
and the variance of k as 



(fc) = 



£+1 ^2_/?+l 



a or 

Substituting the Poisson multiplicity distribution (jHJ and the gamma distribution H15[l into Eq. (|14|l . one finds 

Let us first discuss Eq. Ijltjfl in two limiting cases. When the variance of the source number a 2 vanishes, Eq. i|16|) 
reproduces, as expected, the result i|l(J|) with k — (k). For a 2 — > oo, one gets 

«P)W = o-^ = (P>. (17) 

Thus, the momentum ((p)^)) is independent of A" and it equals the inclusive average - the correlation of the transverse 
momentum and multiplicity entirely disappears. 

The correlation of the transverse momentum and multiplicity is reduced due to the source number fluctuation. 
To quantify the effect of the correlation reduction, Eq. l(TB|) is expanded in the powers of (1 — N/((N)k)), where 
((TV)^) = (k)(N), and only the linear term is kept. Then, Eq. changes into 

(( P )f w) ) = «' + &'(l-^y), (18) 



where a' = a — 1/(N) and 



^ 2 W(fe) 3 + (fe) 4 (19) 

~ (a 2 (N)(k) + {k) 2 -a 2 ) 2 ' 1 ' 



When a 2 0, then b' — » b. For a 2 — > oo, b' — > 0. One also observes that 6' — > 6 for (fc) — > oo. Eq. (|19|l can be 
rewritten using the scaled variance u = a /(k) as 

6', 6. + 1 



(w(^} + l-w/(fc)) 2 
but it does not eliminate an explicit dependence of b' on (k). 



V. MODEL PREDICTIONS 



We discuss here predictions of the formula (|18fl , taking into account the conditions of the NA49 experiment where 
the particles produced in p-p and Pb-Pb collisions are observed, as mentioned in the Introduction, in the forward 
rapidity window [4.0, 5.5]. The coverage in the azimuthal angle is also incomplete. The average multiplicity of charged 
particles in p-p interaction is (N) = 1.4 in such an acceptance |2|. 

The formula (|18() is first interpreted in the Wounded Nucleon Model belonging to the transparency class and, 
as explained in the Introduction, the source number k is identified with the number of wounded nucleons from a 
projectile N pro j. The number iVp ro j is experimentally determined on the event-by-event basis, using the zero-degree 
calorimeter. Due to the finite energy resolution of the calorimeter, the number -/Vp ro j fluctuates with the dispersion 
identified as a. 

In Fig. ^ there is shown the ratio b 1 /b, which measures the correlation reduction, as a function of the average 
number of wounded nucleons (iV pro j) for several values of the dispersion a of the wounded nucleon number. As seen, 
the observable correlation strength is strongly reduced for the peripheral collisions and it does not match to the p-p 
point (b'/b — 1 at (A^j) = 1), as one can naively expect. 

The dispersion of jy n rn j is usually not constant but it changes with (N pTO j). The energy resolution of the NA49 
zero-degree calorimeter [11] which is recalculated into the dispersion of projectile participant number N pm j is 



<W(AW) = 0.77y/A-{N pioj ) - 0.28 - 0.032 (A - (AT proj » , 



(20) 
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FIG. 1: (Color online) The ratio b'/b as a function of the 
average number of projectile wounded nucleons (iVproj} in 
Pb-Pb collisions for several values of the dispersion a of 
Nproj. The most upper line corresponds to <j = 1, the lower 
one to a = 2, etc. 



b'/b 




FIG. 2: (Color online) b'/b vs. (iVproj} in Pb-Pb collisions. 
The upper curve ('res') shows the effect of calorimeter reso- 
lution; the middle one ('mix') represents the mixing model; 
the lower curve ('res+mix') combines the effect of calorime- 
ter resolution and the mixing-model prediction. 



where A is the atomic number of the projectile nucleus. For the lead (A — 208) projectile, the formula l|20l) gives a 
which varies only between 3.5 and 4.0 for 1 < (iVproj) < 150, but it drops to 1.6 for (N pro j) = 200. The prediction of 
the formula l|19l) with the dispersion of the wounded nucleon number given by Eq. I|20l) for A = 208 is represented by 
the most upper curve labeled as 'res' in Fig. [21 

As mentioned in the Introduction, the transparency model, which assumes that the projectile participants mostly 
contribute to the particle production in the forward hemisphere, fails to describe 0,0 the large multiplicity fluctua- 
tions observed in the forward rapidity window |(J . The mixing models 0, Q , where the wounded nucleons from the 
projectile and target equally contribute to the forward rapidity window 1 , agree, at least qualitatively, with the data 
[ElI So, let us consider the model predictions concerning the correlation pt vs. N. The average number of sources 
(k) then equals ((iVp ro j) + (iV targ ))/2 where iV targ is the number of wounded nucleons from the target nucleus. Since 
we consider collisions of identical nuclei (iVtarg) = (Nproj)) and thus (k) = {N pTO j). 

We first neglect the experimental effect of finite calorimeter resolution, and we assume that -/V pro j is precisely 
determined - iV pro j does not fluctuate and (iV pro j) = iV pro j. The number of particle sources, which contribute to the 
forward hemisphere, fluctuates due to the fluctuations of Nt&rg at fixed iV pro j and due to the collision dynamics: a 
wounded nucleon from the projectile (target) is assumed to contribute with probability 1/2 to the forward (backward) 
hemisphere. Then, the variance of the particle sources, which contribute to the forward hemisphere, is a 2 = N pTO - 3 /2 + 
cr 2 arg /4, where <7 2 arg is the variance of -/V targ at fixed N pro j. Using the numerical results obtained within the Hadron- 
String Dynamics model |13| , which were used in 0, , cr 2 arg can be roughly parameterized as 

<U(Nproj) = 3.28iV proj - 0.0158iV 2 roj , 

for the Pb-Pb collisions. As seen, otarg vanishes for -/V pro j = and N pm j — 208, and it reaches maximum equal 13.1 
at iVp ro j = 104. Thus, the variance of the particle sources is 

^(iVproj) = 1-32 iVproj - 0.004iV p 2 roj , (21) 

with maximum equal 10.4 at iV pro j = 165. Substituting (k) — iV pro j and a given by Eq. I|21|l into the formula i|19|) . one 
finds the result shown in Fig. [5] by the middle curve labeled as 'mix'. As seen, the mixing model leads to a sizeable 
reduction of the correlation pt vs. N, and the reduction does not much depened on the collision centrality. 

Finally, we have combined the effect of 'mixing' with that of the finite calorimeter resolution. Treating the two 
sources of fluctuations as independent from each other, the variances have been summed up. The result is shown in 
Fig- Elby the lowest curve labeled as 'res+mix'. 



1 The analysis performed at RHIC en ergi es shows that the contribution of a projectile (target) wound nucleon actually extends far into 
the backward (forward) hemisphere |12| . 
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VI. SUMMARY AND CONCLUSIONS 



The correlation px vs. N, which is evident in p-p collisions, has been studied in the superposition model of nucleus- 
nucleus collisions where produced particles come from k identical and independent sources. At first the number of 
sources has been treated as a fixed number, and then the fluctuations of k have been taken into account. While the 
superposition of particle's sources preserves the correlation strength, the correlation is strongly reduced by the source 
number fluctuations. The calculations are fully analytical. 

The Wounded Nucleon Model has been used to formulate predictions for the NA49 measurement of the correlation 
Pt vs. N in Pb-Pb collisions. Since the measurement, which is underway, is performed in the forward rapidity window, 
it has been first assumed that only wounded nucleons from the projectile contribute there. Taking into account a 
finite energy resolution of the zero-degree calorimeter, which allows one to precisely fix the number of wounded 
nucleons from the projectile, the correlation has been found to be strongly reduced in the peripheral collisions while 
the reduction in the central collisions is quite small. 

Then, the 'mixed' Wounded Nucleon Model 0i@ has been considered. The model, which appears to be successful 
in describing the multiplicity fluctuations at fixed number of the projectile wounded nucleons [H,li|, assumes, contrary 
to a natural expectation [T^| , that the wounded nucleons from the projectile and target equally contribute to the 
forward rapidity window. In this case, the correlation px vs. N seen in the forward rapidity window is strongly 
reduced and the reduction weakly depends on the collision ccntrality. 

The correlation px vs. N seems to be responsible jj| for the transverse momentum fluctuations [2| as quantified by 
the measure 3>(pr) Q The data on the px vs. N correlation combined with the existing data on the multiplicity 
0, |(| and transverse momentum fluctuations - all obtained in the same experimental conditions - will hopefully 
allow one to formulate a coherent picture of the event-by-event fluctuations observed in nucleus-nucleus collisions. 
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